arXiv:chao-dyn/9808010vl 10 Aug 1998 


Boosting Sonoluminescence 


Joachim Holzfuss^, Matthias Riiggeberg^, Robert Mettin^’^ 

^Institut fiir Angewandte Physik, TU Darmstadt, Schloflgartenstr. 7, 64289 Darmstadt, Germany 
^ Drittes Physikalisches Institut, Biirgerstr. 42-44) 37073 Gottingen, Germany 
(Accepted for publication in Phys. Rev. Lett.) 

(Copyright 1998 by The American Physical Society) 

Single bubble sonoluminescence has been experimentally produced through a novel approach of 
optimized sound excitation. A driving consisting of a first and second harmonic with selected 
amplitudes and relative phase results in an increase of light emission compared to sinusoidal driving. 
We achieved a raise of the maximum photo current of up to 300% with the two-mode sound signal. 
Numerical simulations of multimode excitation of a single bubble are compared to this result. 


PACS numbers: 78.60.Mq, 42.65.Re, 43.25.-fy 

By focusing ultrasonic waves of high intensity into a 
liquid, thousands of tiny bubbles appear. This process of 
breakup of the liquid is called acoustic cavitation. The 
bubbles begin to form a fractal structure that is dynam¬ 
ically changing in time. They also emit a loud chaotic 
sound because of their forced nonlinear oscillations in 
the sound held Q. The large mechanical forces on ob¬ 
jects brought into contact with the bubbles enable the 
usage of cavitation in cleaning, particle destruction and 
chemistry. Marinesco and Trillat found that a photo 
plate in water could be fogged by ultrasound. This multi¬ 
bubble sonoluminescence (MBSL) has been analyzed by 
many researchers, and a great amount of knowledge has 
been gained ||]. The discovery by Gaitan Q that it 
is possible to drive a single stable bubble in a regime, 
where it emits light pulses of picosecond duration 
called single bubble sonoluminescence (SBSL), has been 
encouraging scientists to explore the phenomenon and 
the associated effects with a multitude of experiments, 
theories and simulations. The experimental results show 
picosecond synchronicity |^, quasiperiodic and chaotic 
variability of inter-pulse times [^], a black body spec¬ 
trum and mass transport stability [|ll[. The theo¬ 
ries to explain the source of SBSL range from hot spot, 
bremsstrahlung [p^, c ollision induced radiation , and 
corona discharges [|14| to non-classical light El- Numer¬ 
ical simulations have been focusing on the bubble dy¬ 
namics, behavior of the gas content ||^,^, properties in 
magnetic fields Q and the stability of the bubble p^ . 
However, the final answer concerning the nature of SBSL 
still remains open. 

The amount of energy concentration from low energy 
acoustic sound waves to 3 eV photons |||,|l^ raises the 
question, whether the effect can be up-scaled. In this 
paper, we report on experimental enhancement of SBSL 
light production by a bimodal excitation of the bubble 
oscillation. The experiment follows an idea stated in 
[ pOf . We also present numerical simulations of multimode 
sound driving that reveal how multiharmonic excitation 
can adapt to the highly nonlinear bubble oscillation in 


the sense of a strong collapse. 

The experimental setup is as follows: An air bub¬ 
ble is trapped in a water filled cell consisting of two 
piezoceramic cylinders connected via a glass tube [ Elt - 
The levitation cell (“Crum cell”) is standing up¬ 
right with a glass plate covering the lower end of the 
cell. The upper end remains open. A video camera 
pointing from the side allows for online monitoring of 
the experiment. The experiments were done with dis¬ 
tilled and degassed water at room temperature and an 
ambient pressure of 1 atm. The bimodal driving signal 
Pe{t) = Pi cos(27r/t) -I- P2Cos{2Tr2ft cj)) is produced by 
synchronized sine wave generators that allow to fix the 
amplitudes Pi, P 2 , and the relative phase cj). 

Using a multifrequency driving signal however is com¬ 
plicated by two facts. First, the transducers have a com¬ 
plex transfer function; second, the standing wave con¬ 
ditions at each frequency in the cell have to be obeyed 
[p^ . Therefore, multifrequency driving results in space 
dependent phases and amplitude relations and thus in 
an effective sound signal Pa(r,z,t) (with cylinder coor¬ 
dinates r,z of the levitation cell). To measure the am¬ 
plitudes and relative phase that actually appear at the 
bubble position, a small hydrophone is used. The correct 
position is adjusted by hrst focusing the camera on the 
bubble and then inserting the hydrophone at the bubble 
site. The driving signal is digitally recorded and phases 
and amplitudes are recovered via a Fourier transform. 
The light flashes emitted at collapse are measured with 
a photomultiplier. 

Levitating small oscillating bubbles of volume V (t) in 
nonzero gravity is possible through the interaction with 
the driving sound field Pair, z, t) which depends on space 
and time. The time averaged primary Bjerknes force j^ ] 

FB = -{V{t)VPair,z,t)) (1) 

can overcome the buoyancy force and attract the bubble 
to a fixed position in space. Weakly sinusoidally driven 
bubbles of equilibrium radius i?o are trapped near a pres¬ 
sure antinode if they are driven below their linear reso¬ 
nance (Minnaert) frequency /m = (27ri?o)~^\/3 kPo/p ~ 
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3/i?o[Hz] for the experimental conditions used here (with 
polytropic exponent k, ambient pressure poi and liquid 
density p) However, the situation is more compli¬ 

cated for strongly driven bubbles and also for multi¬ 
modal excitation, where the standing wave pattern in the 
resonator, the Bjerknes forces and thus the bubble posi¬ 
tion are changed by a variation of the sound signal pa¬ 
rameters Pi, P 2 , and (j). The bubble oscillation responds 
to the sound signal at the trapping site. 

In the experiment, we proceeded in the following way: 
For fixed drive amplitudes Pi and P 2 , a bubble is in¬ 
jected into the fluid with a syringe. Once the bubble 
fixes itself spatially at a stable position, where the Bjerk¬ 
nes force equals the buoyancy force, the phase difference 
between the locked sine wave generators, one operating 
at / = 23.4 kHz and the other at 2/, is sequentially in¬ 
creased while the sonoluminescence (Sy intensity and 
the bubble itself are monitored. Fig. ^ (lower) shows 
the SL intensity as a function of the phase difference for 
Pi K, 1.25 bars, P 2 « 0.3 bar. With increased phase dif¬ 
ference two maxima appear in the light intensity. The 
dashed line indicates the maximum achievable SL inten¬ 
sity using single-mode driving. This value is obtained 
shortly before and after the two-mode experiment to al¬ 
low direct comparison by keeping all other experimen¬ 
tal conditions unchanged. It is seen that the two-mode 
driving yields 100% more SL intensity than the maxi¬ 
mal single-mode driving. By further increase of Pi and 
P 2 at selected phases an intensity gain of 300% can be 
achieved, as shown by the open circles. Beyond that, the 
bubble gets destroyed. 

Fig. 1^ (upper) reveals, that with increased phase dif¬ 
ference the bubble traverses vertically through stable and 
(surface) unstable regimes. 



Phase [deg] 

FIG. 1. Bubble response for two-mode driving as a func¬ 
tion of the phase difference (in degrees) between the driving 
sinusoidal signal and its second harmonic. Upper: Vertical 
position of the bubble. The thick dotted lines denote un¬ 
stable bubble behaviour. Lower: Photo current. The open 
circles show the maximum SL intensity achieved. The dashed 
line is the maximal photo current for pure sine wave driving. 


Numerical simulations have been carried out using the 
Gilmore model which describes the radial motion of a 
single bubble. The model includes the usual components 
of the Rayleight-Plesset equation like surface tension 
a and liquid viscosity p, and also the compressibility of 
the liquid to allow damping of the bubble motion by the 
shedding of shock waves. 
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R is the bubble radius, Rq = 5 pm its equilibrium radius, 
and C, p, and p are the speed of sound in the liquid, 
its density, and the pressure at the bubble wall, respec¬ 
tively. H is the enthalpy of the liquid. Parameters were 
set to co=1500 m/s, po=998 po=l bar, k=4/3, 

tT=0.0725 N/m, /r=0.001 Ns/m^, n=7, H=3000 bars. 
a=i?o/8.54 is a hard-core van der Waals-term |^. The 
pressure at infinity includes the multimodal driving pres¬ 
sure: Poo = P0+ Pe{t), Pe{t) = Y.m=l PmCOs{2TTmft + 
4^m) • 

First we calculated the driving sound signal that would 
lead to the most violent collapse, indicated by the small¬ 
est minimum radius during a bubble oscillation cycle us¬ 
ing the above equation. The search for suitable pressures 
Pm and phases (pm was carried out by a heuristical opti¬ 
mization algorithm ||^ with the boundary condition of 
a constant driving signal power, i.e., P^ = J2m=i = 
const. Pe was fixed to 1.3 bars and the driving frequency 
was the same as in the experiment. Comparing equal 
power signals is convenient, because the power stays con¬ 
stant upon phase changes, making it possible to compare 
numerics and experiments. 

Fig. 1^ shows the driving pressure and the bubble re¬ 
sponse of different driving signals. A strong increase 
in the maximum radius can be seen already by adding 
just the second harmonic to a sine wave. The numer¬ 
ically computed optimal phase difference is 166.4 deg 
and the individual amplitudes are Pi = 1.026 bars and 
P 2 = 0.798 bar. The radius and the adiabatically cal¬ 
culated temperature around the collapses are shown in 
Fig. ||. It is seen, that the bubble radius at collapse is 
decreased by a large amount and is approaching the van 
der Waals hard core already for the two-mode driving. 
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FIG. 2. Time series of the driving (top) and the radius 
of calculated bubble collapses (bottom) for single (dashed) 
and optimized mnltimode driving signals (two-mode: line, 
eight-mode: dotted). 

Also the maximum temperatures almost double. The 
higher-mode driving signals are better adapted to the 
nonlinear bubble oscillation than the sine signal: They 
show a deeper rarefaction phase before collapse, followed 
by a more rapid rise to the compression phase during 
collapse. 

The calculations for optimal eight-mode driving ex¬ 
hibit only small additional gain compared to bimodal 
driving. Because of the increased difficulties regarding 
the spatial stability of bubbles in the resulting compli¬ 
cated sound held, an eight-mode driving may not be 
worth being considered experimentally. Though two¬ 
mode driving is an early truncation of a series expansion, 
one sees that already this approximation shows a trend 
for a more intense driving of this nonlinear system. 

The optimal results are located on a single broad 
plateau in parameter space. This is in contrast to the 
experimental hnding of two maxima. To understand the 
reason of this obvious discrepancy, the bubble model [Eq. 
(^] has been integrated numerically along with the pri¬ 
mary Bjerknes [Eq. (j^)] and the buoyancy forces to ex¬ 
amine spatial dependencies. This is also motivated by 
the observation that the vertical position of the bubble 
is altered when the phase is changed (Fig. H upper). The 
change of position leads to different effective excitation 
amplitudes for / and 2/ and thus to a more complex sce¬ 
nario. The system of equations is integrated using the 
spatially dependent driving force 

Pe{t) = Pi cos(27r/f) cos (fc/2 z) 

+P 2 cos (27r2/t -I- (j)) cos{2kz — 7r/2) 
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FIG. 3. Zoom into the first collapses of Fig. H Shown 
are the curves for single-mode (dashed), optimized two-mode 
(line) and optimized eight-mode (dotted) driving. The time 
series of the bubble collapses exhibit a decrease in minimum 
radius (left) and increase in the adiabatically calculated tem¬ 
perature (right) as a function of the number of modes in the 
driving sound. The minimum radius comes very close to the 
van der Waals hard core, shown by the horizontal line in the 
left graph. The time axis is shifted so that the collapses take 
place at 0, 0.5, and 1 ns, respectively. 

{k is the acoustic wavenumber 2TTfjco, Pi = 1.25 bars, 
P 2 = 0.357 bar, / = 23.4 kHz). The spatial modes 
are approximately the same as the experimental modes, 
which have been measured with a needle hydrophone. 
The points in vertical z—space where the Bjerknes force 
vanishes and the stability criterion is met represent the 
position of the bubble. The resulting minimum radii are 
shown in Fig. 0 . Comparing this with the experimental 
results in Fig. ^ shows a very close agreement. 





FIG. 4. Numerically calculated vertical bubble position 
(upper) and resulting minimum radius (lower) as a function 
of the phase difference for double harmonic driving of a bub¬ 
ble. The dashed line in the lower plot is the minimum radius 
for single frequency driving with the same power. The solid 
straight line is the van der Waals hard core. 
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The almost sinusoidal variation of the position of the 
bubble gives rise to two minima of the minimal ra- 
dius/phase dependence. Each of these minima is smaller 































than the one of the single-mode driving. The minima co¬ 
incide with the experimental observation of increased SL 
intensity. The slight asymmetry in the experiment can 
be described by the difference in acoustic impedance of 
the glass bottom and the open top of the cell. Also, im¬ 
perfect standing waves may lead to small traveling com¬ 
ponents in the experimental driving. Changing the am¬ 
plitude ratio of the driving signal closer to unity while 
keeping the power constant results in a complex scenario 
of stable bubble positions and effective drivings including 
hysteretic jumps. 

In summary, we have shown that a bimodal sound ex¬ 
citation can enhance light production of SBSL. Though 
spatial modes play a crucial role in double harmonic driv¬ 
ing, it increased the photo current to a gain of maximally 
300% compared to sine excitation. We suppose that mul¬ 
tifrequency driving can shift the bubble oscillation to 
a regime of strong stable SBSL which is not reachable 
by pure harmonic driving. Numerical simulations of an 
acoustically driven bubble including Bjerknes and buoy¬ 
ancy forces show that the increased SBSL light intensity 
is caused by a larger compression. To give quantitative 
estimates, however, elaborate models have to be consid¬ 
ered that include gas dynamic equations for the interior 
of the bubble and thus can model the shedding of a shock 
wave inside a bubble . 

Other methods have been proposed to increase the vi¬ 
olence of bubble collapses. For example, calculations for 
thermonuclear D-D fusion in D 2 O within this context 
have been done using a large pressure pulse superimposed 
on a sine wave Jl^ . However, whether advanced forcing 
by higher modes is large enough to achieve a reasonable 
neutron production rate is an open question. Apart from 
sonoluminescence, the increase of cavitation strength by 
means of optimized multiharmonic sound signals [pot can 
also be of use in the context of sonochemistry |3l[| and 
related areas, where higher reaction rates could be in¬ 
duced. 
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for making the research possible. The work has been 
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the DFG. 
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